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1. INTRODUCTION 

As an important object of analysis and synthesis, you can also specify speech signals that are 
recorded at the outputs of microphones under the acoustic influence of oral speech [1-7]. Their distinguishing 
feature is the non-stationarity due to alternation of speech sounds. Therefore, the processing of finite 
dimension vectors adequately reflects this property. 

The degree of concentration of their Euclidean norms (energies) is often used as the most important 
criterion in synthesis of signals, in some subdomain of the used transformation definition, for example, 
Fourier [8]. Thus, description of the signals properties from standpoint of splittng their transformations 
definitions domains into subregions, is a useful technique for analysis and synthesis of them. Such an 
analysis and synthesis seemto be called subband. 

Subband analysis and synthesis is widely used in framework of the Fourier transform. 
However, one cannot assume that it is optimized from the standpoint of estimating the accuracy of the signals 
subband characteristics. In particular, the existing theoretical foundations of subband analysis /synthesis 
are not sufficiently developed, since only digital filtering issues are sufficiently developed. However, even in 
this case, aspects of informativity are not sufficiently developed from the standpoint of distortions, 
which affect interpretation of the results obtained of the used process ing, for example, accuracy of Euclidean 
signal rate fraction determination in a given subregion (subband) of the Fourier transform definition domain. 
Thus, development of subband analysis/s ynthesis methods from the standpoint of achieving high information 
content from the obtained results is relevant. 
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It is also obvious that the use of Fourier basis, although important, is not the only way to represent 
signals in the corresponding spaces, which make it possible to display the properties of signals that do not 
appear clearly in the time domain. Therefore, it seems appropriate to develop fundamentals of the subband 
analysis and synthesis theory for a certain class of transformations based on orthogonal bases. The solution 
of this problem is the content of the present work. It has been shown that the subband properties of signals 
are adequately described using the subband matrix apparatus, whose eigenvectors are a natural orthonormal 
basis in the space of real vectors of the corresponding dimension. 


2. METHODOLOGY 


In this work, it is assumed that the vector X¥ = (Mies) , where the prime means transposition, 


consists of components obtained either as a result of an equidistant discretization of a continuous signal 
(a function of time). 


x(t), t =[0,T], T = (N — DAt (1) 


In steps At, or from components obtained as a result of synthesis on the basis of a certain criterion, for 
example, during the formation of signal-code constructions in the information transmission systems. 
This refers to compliance, and the boundedness condition ofthe Euclidean norms. 


x, =x(KAt),k =1,..,N (2) 
N 

Ie IEC>) R) <0 6) 
k=l 


The widest used methods of analysis can be attributed to separation of the signal into additive components. 


R 
22 (4) 


Which meet the specified requirements and selection (filtering) of these components. Based on this 
technique, the tasks of extracting components from information signals during multichannel information 
processing [1-3], as well as cleaning signals from the distorting effects of random noise [4, 5], are solved. 
Moreover, note that one of the most important tasks of analyzing empirical data is detection and recognition 
of signals, for which the selection of information components is used as well [6]. 


3. RESULTS 
Fundamentations of the generalized subband analysis and synthesistheorfor symbol F means the set 
in the general case of complex functions, 


F={ ¢,(0), te(ab);n=0,1,... }. (5) 


that satisfy orthonormal conditions, 


Sin = (G49) =| $08, Ow(dt = 6,,. © 


where the asterisk denotes complex conjugation; On, 


is Kronecker symbol; (a,b) is scope, including, 


possibly, the whole numerical axis; w(t) - non-negative weight function 


w(t) = 0,t € (a,b) (7) 
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For simplicity, such functions will be called basic. Examples of such bases commonly used in signal 


processing are the followings: 
Exponential Functions 


Q (t) =exp(-jtn),t e (~m, 7);n =0,... , w(t)=1 (8) 


Cosine basis 


M) =1/ 7,0 (t) =2cos(tn)/ m,t € (—7, 7);n=1,.... w(t) =1 (9) 
Sine base 
Q(t) =0; 0, (0) =2sin(tn)/7,t €(—7, 2);n=1,.... w(t) =1. (10) 


Chebyshevpolynomials the first kind [8, 9] 

& (t) =1/ 2;@,(t) =2cos(n acos(t))/ z,t €(-1,1);n=1,.... wh) =-?)y . aD 
Hermite Polynomials [9] 

¢ (t) =(—1)" exp(t* /2)d exp(—t* /2)/ dt,t € (0, 0);n =0.... 

w(t) =exp(—t / 2) (12) 
Bessel functions [9] 

d (th=J,(u,t/b),t e(0,b);n=1,..; wt) =t / (J7 (u,v 20, (13) 


where 44 ,1=1,...- zeros of the Bessel function the first kind J, V> Q0. In the framework of this work, 


the linear state of the form is called as the transformation of source data (vector (2)) 


X(t) 2 TAO, (14) 


In accordance with the established terminology, the left part of (14) can be called transformant 
of the original data in the chosen basis or simply the transformant. Based on the property (6), 
the representation is valid (inversetransformation). 


x, =| XO¢ wat (15) 


It is also obvious that, based on the property (6), the following equality is held, 


læ I} = w(t)| X(t) P dt. (16) 


a 
Xg = 


N 
k=l 


a w——— S 


Which is a generalization of the Parseval equalty for Fourier transform [10]. It can be represented in 
the following sub-band form. 


Bulletin of Electr Eng & Inf, Vol. 9, No. 3, June 2020 : 964 — 972 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 g 967 
R 
=>, PG); (17) 
r=l 


Where terms are determined by the following integrals 


P(x)= | WOIXOP dt (18) 


tel’. 


For real bases by intervals (subbands) of the form 

T. =(a,,b.),a<a, <b. <b;r=,...,R (19) 
or Symmetric unions for complex bases 

T, =(4,,-4, )UG,»8,) (19) 


It is assumed that the number of such subbands is finite, although their width may be unlimited, for 
example, in the case of using Hermite polynomials. It is obvious that (18) thform integrals determine 
the distribution of Euclidean norms (energies) squares, obtained of the empirical data and can serve as an 
important tool for their analysis. It is not difficult to get their representation directly in domain of the source 
data. To do this, substitute the representation (14) in the definition (18). As a result of simple 
transformations, we obtain the following quadratic form, 


P(E CA (20) 


where C. = {c ],i,k =],.., N - symmetric matrix with real elements 


c= | WOOF Odi 21) 


tel’. 
It is also clear that, since characteristics (18) are not negative, the quadratic forms (20) 


are non-negative, that is, matrices with elements (21) are non-negatively defined. Therefore [11] they can be 
represented as 


C.=QOL,O., (22) 
where Q. is orthogonal matrix of real eigenvectors 

0.0,=Q,0, =I =diagil....,1), (23) 
which is a solution to (24) 

C,Q, =Q,L,. (24) 


But L, = diag (A ,..., Ay )- diagonal matrix of real eigenvalues, which we assume ordered 


in descending order 
A ZA a2 Ok Sla N a (25) 


In the following, matrices with elements of the form (21) are called subband ones. In addition to 
relations (22), (23) and (25), some more specific relations for their own numbers and vectors can also be 
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obtained. In accordance to the definition (24), equality must be maintained for the individual components 


of the eigenvectors Gy. = (Gy son Dry) 
Adk = =$ c Cim Ink i a 


A. Tie = | WA OGL Odt (27) 


tel’. 


Substitution into which representations (21) gives, where the function under the integral is a complex 
conjugate transformant of an eigenvector 


G, (t) = > p Od (28) 


Relation (25) allows to get an idea for projection on the eigenvector (scalar product) of the initial vector 
à ok 
Adi (=A AI=AY ah = | wWOX OG! Oat 29) 
m= teT. 


Note that projection is determined by transformant segments of the eigen and original vectors in the selected 
subband. Similarly, you can get equality 


ZaD AD = | OEO Od (30) 


From here, taking into account the orthonormal eigenvectors (21), we obtain that for the segments 
of the transformants of eigenvectors the folowing equalties are held. Transformants of various 
eigenvectors of 


nk (31) 


have the property of double orthogonality [10] (simultaneous orthogonality throughout the domain and its 
subband), 


| wE (t)GP (t)dt = | w(t)G’ (t)G) (t)dt =0, (32) 


a tel, 


While the eigenvalue determines fraction of square of the norm of the corresponding eigenvector falling into 
the subband 


b 
= | wOIG OP ar<| wOIG,P dt=Ln=l.....N (33) 


tel, 


Note that the equality to the unit of the second integral is obtained on the basis of the Parseval 
equality (17) and the orthonormal nature of the eigenvectors. Bearing in mind the decomposition (22) 
and definition (29), from (20) it is not difficult to get an idea for the characteristic (18) 


N 
PŒ@=) 4d æ, (34) 
k=l 
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which allows us to simplify and in particular, parallelize the procedure for calculating quadratic forms, when 
the eigenvectors are pre-computed for their repeated use. When dividing the original vector into additive 
components (see (4)), you can use the ideal requirement 


YO=XO, t ET; Y O=0,tET. (35) 


where r (t) - transformant vector 


Toa (36) 


Obviously, the requirement (35) can be fulfilled exactly if the subband is equal to the entire domain of the 
transformant. Therefore, it seems appropriate to use the following functionality as a measure for 
the error in its implementation. 


ER,3.)=B | wOIXO-Y.OP dt+d-f) | IXOP dt. 37) 


tel. t¢T, 
where the parameter 


O< B<l (38) 


Serves to set the weight (importance) of the components. You can show the view vector 


y, = B(L- PI +2B-YC,)°Cx. (39) 
In case of equal weights, 
p =0,5, (40) 


relation (39) gives 
y.=CX (41) 


From here, taking into account (21), we obtain the following representation for the components 
of the resulting vector; 


Vir = | WOA OX Odt. (42) 


teT, 


Thus, the components of the resulting vector are completely determined by the segment of the transforma nt 
of the original data in the selected subband. Obviously, this conclusion can be extended to the vector (39), 
since it is completely determined by the vector (41). Based on the definition of (21), we obtain the equality 


y C= (43) 


when subbands are adjacent to each other and cover the entire domain. Obviously, this corresponds to 
obtaining components of the vector that will be analyzed, which meet the requirement of components 
additivity (4). The problem of signal synthesis with specified properties is often considered during 
transmission of information and remote control. The criterion of energy concentration (the square 
of the Euclidean norm) is widely used in a certain subband of the entire frequency band. This criterion is easy 
to generalize in the case of any transformations, and to formulate the problem of finding a vector that 
maximizes the characteristic (20) 


ZC.Z=max X¥C_X, VX € R” (44) 
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The solution of this problem is the matrix corresponding eigenvector with elements (21) corresponding for 
the maximum eigenvalue. 


z=ğ (45) 


It is clear that subband matrices serve as the basis for solving the considered problems of generalized 
subband analysis / synthesis of vectors. In the case of symmetry and non -negative definiteness, they have a 
complete set of eigenvectors that can serve as a basis of space. These bases are determined by 
the type of the transformant and the selected subband of its domain. Therefore, it is natural to call them 
subband. It is mteresting to consider specific examples of subband matrices and their eigenvalues and 
vectors. For some types of transformations, it is possible to obtain analytical relations for the elements of 
subband matrices. 

a. Exponential basis (8). 

The elements of the subband matrix (21) in this casecan be represented in an analytical form. 


c, =2sin(A (i—k))cos(a@,(i—k))/aG—k), (46) 
where 
A. =(b, —a,.)/2;0@, =(b.+a,)/2. (47) 


The main properties of subband matrices and their eigenvectors were considered in [12]. It also 
discusses the features of using this apparatus for solving some problems of analysis and signal synthesis. 
In [13], it was shown that the subband methodology can be successfully applied under conditions 
of complete, in which there is a prion uncertainty in identifying trends of empirical data based on adaptive 
filtering. It does notusea priori postulation for the analytical form of the functional dependence of the trend. 
b. Cosine basis (9) 

The elements of the subband matrix are represented as 


c,, =2sin(A_ (i—k))cos(@.(i—k))/ a —k)+2sin(A, (i +k))cos(a.i+k))/ ai +k) 48) 


Relation (21), after substituting the definition (11) into it and replacing the variable f = cos( z) gives a ratio 


similar to (48), with the difference that instead of (47), the parameters are determined from 
the relations 


A. =(acos(b.)—acos(a,.))/2;@. =(acos(b,) +acos(a,))/2 (49) 


4. DISCUSSION 

Current methodologyfor subband analysis of empirical data and synthesis of signals with desired 
properties is developed regarding the problems needed to be solved from the standpoint of subdividing 
domain of the Fourier transform based on the exponential basis (8) (frequency representations). Frequency 
representations also built developing methods of wavelet analysis quite intensively, for which they use filters 
with finite impulse response, calculated in a quite specific way [14]. This is in favor of the conclusion that 
the sub-band methodology is adequate for the tasks of signal V analysis and synthesis. 

At the same time, for some empirical data and signals, the description using other bases 
(for example Bessel functions) may be adequate when analyzing oscillations of circular membranes [15, 16]. 
Therefore, it is important to develop the mathematical foundations ofsubband analysis and signal synthesis 
from the standpoint of partitioning transformation regions into subbands of other transformations using 
orthonormal bases. This allows for a completely regular way to implement the subband methodology in these 
cases,as well [17-25]. 


5. CONCLUSION 

In this paper, the theoretical foundations of subband analysis are developed from the standpoint 
of dividing the transformation definitions domains (transformants) of the initial sampled data into subbands, 
according to orthonormal bases. It is shown that special matrices, which are naturally called subband ones, 
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are the main elements of subband analysis. It is also indicated that subband matrices are symmetric 
and positive definite. Thus, the set of their eigenvectors is a complete orthonormal basis. It is natural to use it 
as the main tool forsubband analysis / synthesis. A criterion is formulated that defines the optimal subband 
filtering of signals on additive components, and a solution is obtained for the corresponding variational 
problem. It is shown that the selected components have an important property, and they are completely 
determined by segment of the transformant in a given subband. The conditions that allow to restore the 
original signal based on a simple summation of the selected components are revealed. Besides, the relations 
defining the elements of subband matrices for some widely used bases are presented. 
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